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Abstract 

We study the classical limit of the ELPR/FK spin foam models by analyzing the large-distance 



o 

(N 

^J^l asymptotics of the corresponding graviton propagators. This is done by examining the large-spin 

asymptotics of the Hartle-Hawking wavefunction which is peaked around a classical flat spatial 
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. geometry. By using the stationary phase method we determine the wavefunction asymptotics. The 
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CN \ obtained asymptotics does not give the desired large-distance asymptotics for the corresponding 

graviton propagator. However, we show that the ELPR/FK vertex amplitude can be redefined 
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■ such that the corresponding Hartle-Hawking wavefunction gives the desired asymptotics for the 

graviton propagator. 
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I. INTRODUCTION 



Loop quantum gravity is a candidate for a realistic quantum theory of gravity and it rep- 
resents a nonperturbative and background independent way of quantizing general relativity 
However, one of its main problems is finding the classical limit. This is difficult to do 
in the canonical formulation, because there are no appropriate solutions of the Hamiltonian 
constraint. But even if one had such a solution, it would be a complicated expression, and 
showing that its transform to the triad representation has a semiclassical limit which im- 
plies the Einstein equations is a daunting task, see j2|. In the covariant formulation, i.e. the 
spin foam formalism, one can compute the transition amplitudes between the spin network 
states, from which one can infer the spin network wavefunction. However, this will be again 
a complicated expression, and it will be difficult to compute the classical limit. 

In spite of these difficulties Rovelli found a way to study the semiclassical limit indirectly 
[sl . His idea was to consider the graviton propagator within the spin foam formalism and to 
study the semiclassical limit by analyzing the large distance asymptotics of the propagator. 
By using an assumption that the flat-space wavefunction has a specific Gaussian form in the 
spin network basis, Rovelli was able to show that the graviton propagator had the correct 
large distance asymptotics. For more detailed studies and further developments see [443]. 

In jo] it was pointed out that the Gaussian wavefunction which had been used to calculate 
;he graviton propagator asymptotics does not satisfy the Hamiltonian constraint, also see 
8|. Although the physical wavefunction ^1/7(7, Jo) is not a Gaussian, one will still obtain 
the desired propagator asymptotics if "^-yiJjjo) is approximated by the Rovelli's Gaussian 
wavefunction for large spins, i.e. 



^7(i;io) ~ ^(jo)exp 



(1) 



Here 7 denotes the spin network graph, ja is a spin of a link of 7, 9a are arbitrary constants 
and a is a numerical matrix. The parameter jo determines the scale of a triangle area in 
the spin network and can be related to the boundary background metric, see 

However, nobody has investigated whether any viable candidate for the flat-space wave- 
function has the Gaussian asymptotic form ([T]). Note that such an analysis has been recently 
performed in the case of canonical Euclidean loop quantum gravity (LQG) theory It 
was shown that the wavefunction has a Gaussian asymptotics, but it is not of the form 



([T]). This result then imphes that the Euchdean LQG graviton propagator does not have 
the desired large-distance asymptotics. In order to be sure about the implications of this 
result for physics, one needs to perform the same analysis in the Lorentzian case. This can 



be done by using ELPR/FK spin foam models 10|, lUl, since these are the only spin foam 
models that have a Lorentzian formulation and give rise to a LQG theory on the spin foam 
boundary. 

In order to satisfy the Hamiltonian constraint, we will consider a boundary spin network 
wavefunction obtained from the spin foam state sum for a spin foam with a spin network 



boundary. This is a spin foam analog of the Hartle-Hawking wavefunction [12|, and it is 
known that a Hartle-Hawking wavefunction satisfies the Hamiltonian constraint. Guided by 
the construction of the fiat-space wavefunction in the Euclidean LQG case we will intro- 
duce the edge insertions in the boundary spin network in order to simulate the presence of 
the boundary background metric. The large-spin asymptotics of the boundary wavefunction 
will be studied by using the stationary phase method. We will determine the conditions 
necessary for the asymptotics to be of the form ([T]), see Eq. (|20l) . Since the graviton propa- 
gator asymptotics is determined by the jo-dependence of the exponent in ([T]), we will focus 
our attention on the coefficient 5* in f l20|) . which is determined by the Hessian matrix for the 
logarithm of the spin-foam amplitude for the boundary wavefunction. 

The method to determine the jo-dependence of S relies on a nontrivial mathematical 
result formulated in Theorem 1. We obtain that S = 0(1), rather than the desired result 
S = 0(1/ jo), see Eq. ( |27|) . The result S = 0(1) implies that the graviton propagator 
behaves as the distance to the fourth power in the limit of large distances. We will also show 
that the S* = 0(1) asymptotics is a direct consequence of the vertex amplitude asymptotics 
(IBip . which is a common feature of all known spin foam models. 

This paper is organized as follows. In section |TT] we introduce the boundary spin-network 
wavefunction with insertions. In section UTTl we rewrite the wavefunction in the form suitable 
for the asymptotic analysis. Section IIVI is devoted to the analysis of the critical points of 
the wavefunction, which play a major role in the asymptotic analysis. We discuss the 
properties of the stationary point equations and outline a method that can be used to 
solve them. However, it is not necessary to solve explicitly the stationary point equations 
since it is sufficient to use certain properties of the critical points. In section |V] we apply the 
extended stationary phase method to determine the asymptotic behavior of the wavefunction 
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in the large-spin limit. A detailed analysis shows that if certain reasonable assumptions are 
satisfied, the wavefunction will have a Gaussian asymptotics. The width of the Gaussian is 
determined by a complex matrix which is essentially the Schur complement of the Hessian 
of the logarithm of the integrand. It depends in a nontrivial way on the scaling parameter 
jo- In order to be able to compare the wavefunction asymptotics with the Gaussian from 
([1]), we need to determine the scaling of the Schur complement in the limit jo — ?■ oo, which 
is done in section I VII An explicit calculation of the Schur complement will not be possible, 
but it will be possible to determine its scaling dependence on jo- Surprisingly, one finds 
that in the leading order the Schur complement scales as a constant in the limit jo — ?■ oo, 
in contrast to the assumed 1/jo scaling in ([T]). This implies that the corresponding graviton 
propagator does not have the distance scaling corresponding to a graviton propagator from 
general relativity. In the final section lylTl we discuss the possible ways to solve this problem 
and to recover the desired scaling of the propagator. It turns out that the most promising 
method is to redefine the vertex amplitude of the spin foam model, and we propose two 
ways to do that. The appendices lAl [B], [Cl |D] and |E] contain derivations of the results that 
were used in the main text. 



II. THE BOUNDARY WAVEFUNCTION 



A boundary state for an ELPR/FK spin foam model can be constructed in the fol- 
lowing way. We expand in the spin network basis I7, j;, Lp), where 7 is the boundary spin 
network graph, j/ are the spins of the edges of 7 and tp are the corresponding intertwiners. 
We then expand each I7, j;, Lp) in the coherent state basis \'y,ji,npi), see [iSj, so that 

1^) = XX / n^^^P' '^-iUunpuio) \-f,ji,npi) . 
7 ji {pi) 

The coefficients ^{j,n\ Jq) are constructed as boundary spin-network wavefunctions with 
edge insertions. The edge insertions are introduced in order to provide the wavefunction 
with the information about the boundary background metric. This is done through the 
background spin parameter jo, so that 



'^^iiunpuio) = \{^^l{JuJo)df{]l) X / \{d''nef\[df{kf)\[W,{kf 



nef,Ji,npi) 



(2) 



dcr=J 
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where a is a 2-complex whose boundary one-complex is 7. The face labels k and the edge- 
face labels n of the corresponding spin foam are fixed to be ji,npi at the boundary spin 
network, kf is a non-boundary spin, which labels a face /, while a unit vector labels 
an edge e and the face (e/) adjacent to e in the 2-complex a. One can also include a sum 
over various 2-complexes a that have the fixed boundary 7 in ([2]), thereby implementing the 
"sum over triangulations" idea. However, this will not affect our analysis, so that we will 
work with a single a for a given 7. The expressions for the face and the vertex amplitudes 
df and can be found in Q, O, 14, 15| and we do not write them explicitly because we 
will need only their asymptotic form for large spins. 

The motivation for the introduction of the edge insertions jo) comes from the con- 
struction of the flat-space wavefunction in the Euclidean canonical LQG [2]. This wavefunc- 
tion solves the Hamiltonian constraint and it is given by a state-sum similar to (|2]). The 
parameter Jq is proportional to the areas of the triangles determined by the background 
geometry triads. The boundary spin network insertions are arbitrary functions of the edge 
spins and jo- Therefore the expression ([2]) is a natural generalization of the Euclidean wave- 
function from {2! to the Lorentzian geometry case. Furthermore, since ([2]) is constructed as 
a Hartle-Hawking wavefunction for a boundary spin network for the ELPR/FK spin foam 
model, it will satisfy automatically the corresponding Hamiltonian constraint and the in- 
sertions will insure that it is peaked around a flat spatial geometry. Since the insertions 
fJ'i{ji,jo) can be arbitrary functions, one can try to choose them such that the asymptotics 
(II]) is obtained. 

The vectors fief are in general defined up to arbitrary phase factors. These phase factors 



can be chosen such that they insure nice g. 
triangulation dual to cr. As discussed in 



uing properties of neighboring simplices in the 
16| . such a choice will fix the phase factors on the 
spin foam boundary 7, and thus give rise to the phase term in ([1]). However, these phase 
factors disappear when the graviton propagator is calculated in the standard canonical 



formalism, see 



61], and therefore their values will not be important for our purposes. 



III. ASYMPTOTIC ANALYSIS 

The wavefunction ([2]) does not necessarilly have the large-spin asymptotic form ([T]). In 
what follows, we are going to study its large-spin asymptotics, in order to find out is there 
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a choice of the insertions such that the asymptotics ([T]) is obtained. If the asymptotics 
is indeed of the form ([1]), the wavefunction ([2]) can be a good candidate for a flat-space 
wavef unction. 

We begin the analysis of the large-spin asymptotics of ([2]) by deflning the large-spin limit. 
Namely, we are interested in the limit 

3 1 = joji, jo oo. (3) 

Here J; G No/2 are spins which are flxed, while jo is the large parameter. 

It is important to note that the scaling of boundary spins ji via the parameter jo will 
induce a similar scaling in some of the internal spins kf, due to the triangle inequalities built 
in the vertex amplitude W^- However, not all internal spins need to be scaled, depending 
on the combinatorics of the two-complex a. The domain of summation in ([2]) will contain 
sectors where all spins are scaled and sectors where only some of them are scaled. Those 
internal spins which must scale do so by a prescription analogous to ([3]). 

The flrst step in flnding the asymptotic behavior of ([2]) is to approximate the sums over 
the internal spins kf with integrals. The wavefunction (|2]) can be then approximated as 

^tU'^p'! jo) ~ A(j^^p«; jo) = 

= f Yldkf f JJ rf^^e/ e^"^^^''^'"'^") , (4) 
Jd J. J 

where the function F is given by 

F{j,k,n;jo) = — log/ii(j/; jo)rf,(ji) + 
Jo V 

+ 1 5^ log dfikf) + - 5^ log W,{j, k, n) . (5) 

D is the domain of integration over spins k and the form (jl]) is suitable for the stationary 
phase approximation. Note that the vertex amplitude Wy is complex-valued in general, so 
that the logarithm is deflned up to a multiple of 27ri. However, this constant factor does not 
influence the subsequent analysis and we can ignore it. Also note that the insertion functions 
Hi depend exphcitly on jo, while df and may depend on jo only through boundary spins 
j and those internal spins k that are constrained to scale via triangle inequalities. 
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We will use the extended stationary phase method [171] in order to approximate the 
integral (jl]). The method will be applicable if the function ([5]) satisfies 

F{j,k,n;jo) = 0{l), (6) 

for jo — oo- This condition will be satisfied on a subset of D where the asymptotic formulae 



for the ELPR/FK vertex amplitude W^, derived in [16|, |20|, l2l|, are valid. See the appendix 
[B]for the explicit expressions. 

When the boundary spins ji are large, i.e. ji = 0{jo), then the integration domain D 
will contain spin foams whose spins are all large. D will also contain spin foams where some 
of the spins are large and other are small. This structure is a consequence of the triangular 
inequalities among the spins which form a spin-foam vertex (rules for the addition of angular 
momenta). Let -Dndg be the set of spin foams in D such that each spin foam from -Dndg 
contains at least one vertex with all spins large. Then D^g = D \ -Dndg is the set of spin 
foams where every vertex in a spin foam from D^g contains a small spin. Consequently 

T — rndg _|_ rdg 
-'7 — -'^ 7 ' 

where J"'^^ and /^^ are defined by taking the integral (jlj) over the domains -Dndg and -Ddg, 
respectively. 

It is not known whether the function F satisfies the condition ([6]) on -Djg, since the 
asymptotic formulae for Wy when some of the vertex spins are large and the other are small 
are not known. On the other hand, the asymptotic formula for in the case when all the 
vertex spins are large is known, see (IBip and (lB2p . so that it can be shown that F satisfies 
the condition ([6]) on -Dndg- This is true because every spin foam from -Dndg contains at least 
one vertex with nondegenerate asymptotics, and therefore the contribution of such a vertex 
to F is given by 

- log W:'^ ^ - log f iVf e^"^^"' + iV(")e-^°^«^) +o(^-^] . 

Jo Jo V J \ Jo J 

Since N±. ^ 0, then 

- log Wf^ ^ + 1 log f ivf + iv(")e-2-^«') + O f ^1 . (7) 

Jo Jo Jo ^ ^ V Jo / 

According to (1A2I) the Regge action S'^^ is of 0(jo), so that the first term in ((Tj) is of 0(1). 

Since the coefficients A^'^."-' are of 0(1), the second term in ([7j) is of 0(j(7^) . Therefore, a 

nondegenerate vertex gives an 0(1) contribution to the function ([5]). 



A degenerate vertex from -Dndg can give a contribution to F of 0(1) or lower, depending 
on the type of degeneracy of each particular vertex. The sum over the insertion functions /x/ 
in F can be chosen such that it is of 0(1). One particularly useful choice for the insertion 
functions is 



f^iUhJo) = exp 



(8) 



Jo 

This choice is very natural for our purposes, since it enforces the flat background metric in 
the boundary state and it gives an 0(1) contribution to F. As far as the the sum over the 
face amplitudes in F is concerned, it is of O(jo^^lnjo), which is subleading to 0(1). This 



is because df{j) is of 0(jo), where g = 1 or g = 2, see [l5| for a discussion of the various 
proposals for df{j). Therefore F = 0(1) on -Dndg? provided that there is no cancellation of 
0(1) terms. Hence one can use the stationary phase approximation for the integral J'^'^s. 

As far as the order of F on D^g is concerned, it can be of 0(1) if the choice ([8]) is used. 
However, the stationary phase approximation cannot be made because of the absence of the 
asymptotic formulas for the degenerate vertices. 

Also note that the extended stationary phase method is directly applicable only if F is a 
Morse function, which means that its Hessian matrix does not have zero eigenvalues at the 
critical points. However, in our case the Hessian of F may happen to be degenerate, so that 
we need to take this fact into account when applying the stationary phase method. This 
will be discussed in detail in section |Vl 

IV. CRITICAL POINTS 

The idea of the stationary phase method is to approximate the integrand in the nonde- 
generate piece of as a sum of Gaussian functions, where each Gaussian is centered around 
a stationary point {j*,k*,'n*) of e^"^. As jo ~^ oo, only the immediate neighborhoods of 



the stationary points will contribute to the integral 17|. Furthermore, only the stationary 
points for which 

ReF(j*,F,r;jo)=0, (9) 

will give a noticeable contribution. The stationary points which satisfy (Q are called the 
critical points. 

Note that the stationary points of e^'^^ are the same as the stationary points of F. 
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Therefore, the stationary point equations are given by 

dF dF dF 

ai-'- wr'' ^r'- J'"' 

The geometric interpretation of these equations has been studied extensively in [18| for 
the Euclidean theory and in [iqI for both Euclidean and Lorentzian versions of the theory. 
In this paper we do not need to go into the details, just let us mention that the condition ([9]) 
and the the stationary point equation for fief are satisfied only if one considers a spin foam 
which is dual to the triangulation of a Regge geometry. As far as the j and k equations are 
concerned, we can unify them by using a common label Xa for ji and kj. Then ^ = gives 

v&dg I f 

where we have used (]A3p . When Xa = kf then the /i; terms are absent in (ITTl) . while for 
= ji the last sum in f lTTj) goes only over the boundary /'s. 

The jo-dependence of (fTTl) is the following. The first sum in (fTTl) is of 0(1), and it gives 
the dominant contribution to the equation. The second sum in flTTl) is of 0(1/ Jq) since the 
factors A^^"'' are of 0(1) and consequently their derivatives with respect to j and k are of 
O(l/jo)- The fourth and the fifth sum in (fTTj) are of 0(1/ jo), while the third sum is at most 
of 0(1). Given this, the critical point equations @ and ffTOj) can be solved by restricting to 
a Regge triangulation and writing a stationary point x* as 

X: = CaJo + da + 0{l/jo), (12) 

where Ca and da are coefficients to be determined. The equations can be then expanded into 
apower series in 1/jo and solved order by order for c and d. This has been done explicitly in 
91] for the case of Euclidean LQG fiat-space wavefunction. In that case the vertex amplitude 
is the 6j symbol, and some explicit solutions can be found. 

However, for the purposes of this paper, it is not necessary to construct explicit solutions 
of (Q and (fTOj) . Rather, we only need to assume that they have at least one nontrivial 
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solution (j*, k*,n*) which is a critical point. If there are no such solutions, the integral I^'^^ 
will be of oII/Jq) for all n > 0, and thus it will not have the asymptotic form ([1]). 

Note that Q^^ = for all / is a leading-order solution of ( TTT]) . since if we neglect 0(1/ jo) 
terms we obtain 

However, such solutions have to be discarded because the corresponding {j*,k*,n*) do not 
satisfy the triangular inequalities. The reason is that the angles 6^^^ are exterior dihedral 
angles of a 4-simplex dual to v. Since a 4-simplex is a convex body, its exterior dihedral 
angles cannot all be equal to zero. 



V. EXTENDED STATIONARY PHASE METHOD 

We are going to determine the large-spin asymptotics of J"*^^ by using the extended 
stationary phase method. As explained in the previous section, we will assume that there 
is a dominant critical point and such a point must satisfy 

e?) ^ , (13) 

for some / and all v. 

As we have already pointed out, the function F may not be a Morse function, and 
consequently we cannot apply directly the well known results \v\ so we will perform the 
calculation step by step. 

Let J = {jiyUpi), K = {kfjfief) and x = {ji, kj, flpi, fief)- We first approximate the 
integrand e^°^ with a sum of Gaussian functions, each centered around a critical point x*. 
The corresponding exponents are obtained by expanding F into a power series around each 
X* up to quadratic terms. The integral I^'^^ then becomes 

X* 

= 5]e^»^*r(J,x*). (14) 

X* 

Here F* = F{x*) = ilmF{x*) is the evaluation of F at the critical point x*, the sum goes 
over the set of all distinct critical points, and A is the Hessian matrix of joF evaluated at 
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X 



A 



ab = JO- 



(15) 



dxadxb 

The scaling parameter jo has been absorbed into A for convenience. 

In order to perform the Gaussian integrations in (1141) . we will split the A matrix into J J 
JK and KK blocks, which will be denoted as A, N and M, respectively 



A N 
M 

Let us rewrite the exponent in a Gaussian integral as 



(16) 



ix — X*) Aix — X* 



1 



{J- rfA{j -J* 



+ -{K - K*y M{K - K*) + {J - J*y N{K - K*) . 

The first term is independent of K and can be moved in front of the integral. By making a 
change of variables K = QK and by making a suitable choice of the matrix Q, the matrix 
Q'^MQ becomes diagonal. Then 



. J Da 



where rua are the eigenvalues of the matrix M, na = [{J — J*)'^NQ'j^, and Da is the one- 
dimensional domain of integration, determined by D and the change of variables K = QK. 
Let us now discuss the integral I*. 

• If nia 7^ and Rem^ < 0, the integral converges. Since = 0{jo), in the limit 
jo oo the result is independent of the domain Da in the leading order, and can be 
written as 

L 



-A. 



27r 



-rria 



1 + 



1 



Jo 



Note that quadratic dependence on generates a term of type (J — J*)^ in the 
exponent, which gives us the desired Gaussian asymptotics. 

If nia 7^ and Rem^ > 0, the integral diverges exponentially in the limit jo ^ oo, so 
that in this case the Gaussian asymptotics cannot be obtained. 

If rria = and Ua 7^ 0, the integral might or might not converge, depending on whether 
the domain Da is compact or not. However, even when it converges, the result will be 
a non-Gaussian function of J — J*. 
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• If rria = and Ua = 0, the integral converges if Da is compact. Most importantly, in 
this case the result is independent of J — J*. Namely, if we denote Da = [aa,Pa], we 
have 

-^a Pa (^a = -^a ■ 

The integrals of this type do not influence the propagator asymptotics. 

Therefore the integral /* will be a Gaussian function of J — J*, if the following conditions 
are satisfied: 

• all nonzero eigenvalues of M have their real part negative or zero, 

• the matrix is projected to zero on the kernel of M, and 

• the domain of integration over the kernel space of M is compact. 

These conditions imply that we can always make a change of variables such that the 
matrices M and are given by 

""m 





M 



N 



N 



where the matrix M is invertible and has a negative-definite real part, and the horizontal 
dimension of the zero-block in N is equal to the corresponding dimension of the zero-block 
in M. 

We note here that these zero blocks appear because the function F has continuous sym- 
metries, which give rise to manifolds of stationary points instead of discrete sets of stationary 
points. As we have shown, integrating over these manifolds does not affect the propagator 
asymptotics, as long as they are compact (otherwise the integral will diverge). A similar 
situation was encountered when applying the extended stationary phase method to the case 
of an Euclidean spin foam model ISj, as well as to the case of a spin foam consisting of a 



single vertex la, l2Q|, l21 |. 



Consequently we obtain 

X* 

where 



l(J-J*)Ts(x*,jo)iJ-J*) 



(17) 



-2nY 



detM 



1 + 



Jo 
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r is the dimension of M , while the product is taken over a for which = 0. Also, 

S = A- NM-'^N^ (18) 



is the Schur complement [22| of the (regular minor of the) Hessian matrix A. 

Note that the asymptotic form f|T7|) is a sum of many Gaussian functions, while the 
desired asymptotics ([1]) is just a single Gaussian function. The expression ( IT7|) can yield a 
single Gaussian if there is a dominant critical point Xq such that |^(xo)| dominates any other 
|^(a;*)| when jo oo. This point can be determined as the one for which the dimension r 
of M is minimal, i.e. when M is maximally degenerate. Consequently 

when jo oo. Finally, the quadratic form in the exponent of ( !T9l) can be decomposed into 
a sum of jj, jn and nn terms. Since j = 0(jo) and n = 0(1) then the jj terms will be 
dominant in the limit jo ^ oo. Therefore 

lf^{r,Jo) ^ A{jo) e^(^-^*)"^(^o)(^-^*) , (20) 

where S is the jj block of the Schur matrix S. Note that ( [T2|) implies j* = 0{jo), so that 
( EUj) can be written as 



a,b 



(21) 



/f*^(j;jo)^^(jo)exp 
From ([20D it follows that 

^ ^(jo) e5(^-^*)^^(^°)(^-^*) + . (22) 

In order to obtain a single Gaussian asymptotics we need to assume that J^^ has a sublead- 
ing asymptotics to that of 1^*^^. Since we do not know how to calculate the asymptotics of 
J^s, there is a possibility that /^^ has the right asymptotics which is dominant with respect 
to However, we will argue that such a possibility would give a classical limit whose 

spacetime geometry has the curvature which varies greatly on small scales, while the corre- 
sponding propagator is that for a flat spacetime, see section IVIII Therefore we will assume 
that implies 

^ A{jo) e^(^-^*)^^(J")(J-J*) . (23) 

However, the asymptotics (l23l) is still not the desired asymptotics ([2]). We need to to 
determine whether or not S = 0(1/ jo). We will analyze this problem in the next section. 
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VI. CALCULATION OF THE EXPONENT FACTOR 



The asymptotic form (!23|) will give the desired asymptotics if S* = O(l/jo)- Note that it 
is very difficult to calculate the matrix S explicitly. However, we only need to calculate the 
leading jo-order of S. This can be done by using the following theorem 

Theorem 1. // the matrix S is nonzero, and if the leading order contribution to A comes 
from terms, we have 

S = 0(A) 

for jo —7- oo. If the matrix S is zero, the wavefunction asymptotics is non-Gaussian. 

The proof is essentially based on the Schur determinant formula, det A = det S det M, 
see [2^, and is given in Appendix O (see also [9]). 



The asymptotic dependence of A on jo can be determined quite easily, if the large spin 
asymptotics of the vertex amplitude Wy is known. From (|5]) and f|T5l) it follows that 

^ dxadxb ' ^^^^ 

where V G {/, /, f}, Ai = ^i, Af = df, A^ = W^, and the derivatives are evaluated at the 
critical point x*. Each term in ( HM contributes to the asymptotics of A with some power 
of jo, so that the leading order asymptotics of A will be determined by the highest power 
of jo- 

The insertion functions can be chosen arbitrarily and therefore can give any desired 
contribution of 0(jo). However, they only contribute to the diagonal elements of A, since 
each insertion function fii depends only on the spin of its link, ji. For the choice one 
easily gets from (IMl) 

d'^hgm 25abSal q(}_ 



dxadxb jo \jo. 

The face amplitude df is commonly chosen to be df{j) = 2jf + 1, see [15|. Substituting 

into (|24|) . we obtain 

dHogdf 46abSaf ^fl_ 



dxadxb {2xf + iy Vjo. 
and this is also a contribution to the diagonal elements of A. Note that other choices for df{j) 

have also been proposed in the literature, see for example However, all the proposed 

choices satisfy df = 0{j'f), where g > 1, so that one obtains an 0(jo^^) contribution. 
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Finally, the main nontrivial contribution comes from the vertex amplitude W^. The 
asymptotics of the degenerate configurations of the vertex amplitude is unknown, and such 
vertices can in general give a contribution to A of order 0(1) or smaller. However, the 
asymptotics of the nondegenerate vertices is well studied, see Appendix [Bl Furthermore, 
each spin foam in -Dndg contains at least one nondegenerate vertex. By a straightforward 
calculation one obtains from (IMI) and (]B2p 

log 



vab 



(2) 
vab 1 



(25) 



dxadxi, 

where the three terms on the right-hand side represent contributions of order 0(1), 0(jo^^) 
and 0(j(7^). In the non-Regge cases the contributions are of 0(l/jQ) or subleading. The 
leading term is 



JO) 
vab 



« Oa^b&v^a ^b 



(26) 



while explicit expressions for A^^^^^ and A^^^^ are given in Appendix [El 

The equation ( l25l) is evaluated at some critical point x* . As we have discussed in Sec- 
tion IIVI there are no critical points where all the angles Oi^'' are zero. Moreover, for the 
ELPR/FK vertex amplitude the coefficients n'^'^ and N^'' are also always different from 
zero, so the square bracket in fl26l) is also nonzero. Therefore we see that the term fl26l) is 
always nonzero. Hence the dominant contribution to the Hessian A comes from the vertex 
amplitudes whose spins form a Regge geometry, and it is of order 0(1). Consequently, the 
assumptions of Theorem 1 are satisfied, so that 



5 = 0(1). 



(27) 



The implication of (1271) for the large-distance asymptotics of the graviton propagator 
can be seen from the following result. Consider a generalized Rovelli asymptotics for the 
wavefunction 



^7(i;io) ~ ^(jo)exp 



CfeJoJ 



a, 6 



(28) 



where p > 0. Note that the obtained result (p71) corresponds to p = while the Rovelli 
ansatz corresponds to p = 1. The propagator asymptotic scaling with spacetime distance 
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\x — y\ can be determined by repeating the calculation done in [3|, also see |6[. Therefore 
one obtains for large distances 



^/ N const , , 

(29) 



where G denotes the diagonal components of the graviton propagator. 

The equation (!29l) gives for p = 1 the propagator asymptotics consistent with general 
relativity, while for p = it gives 



G{x,y)^j^^^. (30) 



The asymptotics fl5U|) is not consistent with general relativity. 

VII. DISCUSSION AND CONCLUSIONS 

The result ( 1271) has been derived under certain assumptions, so that one would like to 
know is it possible to relax the assumptions such that the desired classical limit is obtained. 
The first thing one can try is to change the insertion functions fii, since these functions can 
be chosen freely. The insertion functions could be chosen such that they cancel the 0(1) 
terms in the Hessian A. However, these functions can only change the diagonal elements 
of A, while the off-diagonal elements will still have the 0(1) terms. Note that we have 
introduced the insertion functions in the simplest possible way, namely as multiplicative 
factors for the amplitude of each link on the boundary spin network. In the most general 
case a fii can be a matrix function, see 2(], so that this gives an additional possibility to 
change the 0(1) behavior. This possibility should be explored, but the problem is that it is 
difficult to analyze. 

Note that we have assumed that the dominant contribution to the asymptotics of the 
wavefunction comes from a non-degenerate spin foam. The reason was that only in that case 
we know how to calculate the asymptotics. Hence there is a possibility that the dominant 
contribution comes from a degenerate spin foam and that this contribution is such that it 
gives the desired propagator asymptotics. However, there is a problem with this. Namely, 
a degenerate spin foam is such that its every vertex has at least one small spin. This means 
that the corresponding spacetime geometry has the curvature which grately varies at small 
scales, which is not consistent with the propagator asymptotics for a fiat spacetime. 
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The only remaining possibility is to modify the ELPR/FK vertex amplitude W{j,n). 
Note that the 0(1) contribution to A is given by fl26l) . and it vanishes if one of the coefficients 
is zero. Consequently, if the modified vertex amplitude n) had the asymptotic 
behavior 

WU,n)^^^, (31) 
where V{j) is the function from flB2p . then it is easy to show from (^M that 

S = 0(1/ Jo) . 

By using ([8]) for the insertion functions, one would then obtain the correct graviton propa- 
gator asymptotics. Note that W gives a state sum which for large spins looks like a path 
integral for Regge discretization of general relativity, because W has the asymptotics (I3T1) . 
This explains why W gives a graviton propagator with a good asymptotics. On the other 
hand, the presence of the complex conjugate term e~*"'^« in flB2|) gives an unnatural path 
integral, so that it is not a surprise that the corresponding propagator has wrong asymp- 
totics. 

Note that all known spin foam models have the vertex amplitude asymptotics which is a 



linear combination of e^''*"'^^' terms, see [l^ for the Euclidean ELPR/FK model or 23l. l24 1 
for the Barret-Crane model. Consequently one will obtain 5* = 0(1) for the large-spin 
asymptotics of the boundary wavefunction, because the calculation is the same as the one 
presented in this paper. It is also instructive to compare our result with the results of the 
similar calculation for the Euclidean theory done in 18|]. This is done in Appendix [Fl and 
supports our result fl271) . 

There are two ways to interpret the result (1271) . One way is to say that the choice (E]) for 
the boundary wavefunction is not the most general one, and there may exist another solution 
which could give the correct asymptotics. This is of course a possibility, and it is an open 
problem for future research. However, such an interpretation of our result essentially brings 
us back to the original problem of finding a wavefunction which satisfies the Hamiltonian 
constraint and has the asymptotics ([T]). However, it is difficult to see what would be an 
alternative construction to the one we used. 

The other possibility is to use the same construction for the wavefunction and to modify 
the ELPR/FK vertex amplitude, so that the result ( !27|) is circumvented. As discussed 
above, the way to achieve this is to construct a new vertex amplitude which would have the 
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asymptotics of the type (13T1) . For example, if 7^ then one can define the new vertex 
amphtude W{j, n) as 



^= \r2_^2 . (32) 



N+W -N^W* 

where W* is the complex-conjugate of the ELPR/FK vertex amplitude W. The new ampli- 
tude will have the asymptotics fl3T]) . A more general redefinition, valid for A^_|_ = case, 
is given by 

--^(--y^^^^)- (33) 

This expression also gives the asymptotics fl3T|) . Hence the spin foam model defined by the 
new amplitude W will give the correct propagator asymptotics and it will represent a good 
candidate for a spin foam model whose classical limit is general relativity. 

Note that the correct asymptotics of the graviton propagator does not guarantee that the 
classical limit of a spin foam model is general relativity. Namely, the graviton propagator for 
a boundary state is defined as a 2-point correlation function. However, in order to determine 
the corresponding semiclassical equations of motion one needs the effective action, which 
is the generating functional for all n-point correlation functions. Knowing just the 2-point 
correlation function is not sufficient, so that one needs to compute the effective action and 
to show that its classical limit is the Einstein-Hilbert action. 
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Appendix A: The Regge action for a 4-simplex 

The Lorentzian Regge action for a 4-simplex dual to vertex v is given as 

S^^\k) = J2kM;\k) (Al) 

Here kj are 10 spins labeling the faces, while each <d^^\k) is the exterior dihedral angle 
between two tetrahedra of the simplex dual to v which share the triangle dual to /. 
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If all spins kf are uniformly scaled as kf = jokf, in the limit jo ^ oo the Regge action 
scales as 

S^^\k) = 0{jo), (A2) 

since kf = 0{jo) and ef\k) = 0(1). 

Also, if we take the derivative of the Regge action with respect to some spin ka, we obtain 



The first sum reduces to gI^^ if a G f , and is zero otherwise. The second sum is identically 
zero due to the Schlafli identity, so we have 

^ = ^-^©i''^ • (A3) 

Note that this derivative scales as 0(1) in the limit jo oo. Also note that for the 
nondegenerate 4-simplex all dihedral angles G^''' are different from zero, since the 4-simplex 
is always convex. These properties are essential for the derivation of our results. 



Appendix B: Asymptotics of the ELPR/FK vertex amplitude 



The asymptotic properties of the ELPR/FK vertex amplitude were investigated in 



depth in 
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2l|, and neatly summarized in 25|. 



A single vertex amplitude Wy is a function of 10 spins kf and 20 normals Ugf. Some of 
the spins may be scaled as kf = jokf, while others do not scale. In the limit jo — )■ oo, the 
asymptotic behavior of can be split into several cases, based on the possible choices of 
these variables. These are 

1. The nondegenerate case 

In this case we assume all 10 spins scale with jo, and the boundary of the corresponding 
4-simplex has the Regge-like geometry. In the case of the Lorentzian version of the 
theory, the vertex amplitude has the asymptotic formula 



Jo 



(Bl) 



Here a = 1 for the 4-simplex with a Euclidean geometry on the boundary, while a = 7 
in the Lorentzian boundary case, where 7 is the Immirzi parameter. The constants 
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are different from zero and S^^\k) is the Euclidean/Lorentzian Regge action 

m- 

2. The degenerate cases of zero 4-volume 

These are the cases when all 10 spins scale with jo, but the boundary of the 4-simplex 



21| where 



does not have Regge-like geometry. The vertex asymptotics was analyzed in 
it was determined that it has the form 

y y V .19 

Jo 

if the boundary is a 3D vector geometry, while 

= oUo''), WK > 0, 

in all other situations with zero 4-volume. All these cases contribute with zero measure 
in in the integral (jl]) and can be ignored. 

3. The degenerate cases of non-zero 4-volume 

These are the cases when only some of the 10 spins scale with Jq, while others are kept 
fixed. These situations have not been analyzed so far, and the vertex asymptotics in 
these cases is still unknown. Note that such configurations contribute with non-zero 
measure in the integral (jlj), and thus cannot be ignored. 

It is important to emphasize that explicit dependence of the asymptotic formula on 
normals n^f is lost in f lBip . and the asymptotic expression on the right-hand side of f lBip 
depends only on 10 spins kf. Namely, the assumption of Regge-like geometry of the 4- 
simplex implies that its triangle areas kf and its normals fief are fully determined by its 
10 edge lengths /j, which also induce the Lorentzian/Euclidean signature of the metric in 
the 4-simplex. However, given that the number of triangles in a 4-simplex is equal to its 
number of edges, the functions kf{li) can in a generic situation be inverted, and edge lengths 
regarded as functions of the triangle areas. This is possible always except in some particular 
cases where the Jacobian of the transformation is singular. Nevertheless, these singular 
cases contribute with zero measure in the integral (jlj) and can thus be ignored. Given the 
inverted functions h{kf), one can also express the normals Uefik) as functions of kj, which 



therefore remain the only independent variables in (IBip . 
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Note that the asymptotics (IBip can be rewritten as 



WU,n) 



1 



(B2) 



for j — !■ oo where ^(j) = 0{j^'^). 

Appendix C: Proof of Theorem 1 

Here we give a proof of Theorem 1 used in the main text. Let us repeat the statement of 
the theorem, for completeness. 

Theorem 1. If the matrix S is nonzero, and if the leading order contribution to A comes 
from terms, we have 



for jo —7- oo. If the matrix S is zero, the wavefunction asymptotics is non- Gaussian. 

The proof goes as follows. Begin by noting that the Hessian matrix A is non-diagonal. 
Namely, looking at its definition (1151) and the action iQ, we see that the insertion functions 
fii and face amplitudes df contribute only to diagonal terms in A, since each of them is a 
function of a single variable. In contrast to this, the vertex amplitudes Wy are functions of 
10 or 30 variables each, according to the combinatorics of the spin foam 2-complex and the 
possible degeneracy of W^. Therefore, each vertex amplitude will contribute to the diagonal 
terms of A, and in addition also to some off-diagonal terms on each side of the main diagonal, 
in such a way that in every row and column there will be some nonzero nondiagonal elements 
present. 

We want to discuss the dependence of det A on jo in the limit jo oo. For simplicity, 
in what follows we shall assume that det A 7^ 0, and we shall discuss the singular case later. 
The determinant of A is by definition given as 



where p is the permutation of indices 1 . . . R, and R is the rank (and simultaneously the 
dimension) of A. In this sum, there will be some terms which contain diagonal terms of 
A, and terms which do not contain any diagonal element. The first set of terms will have 
contributions of fii, dj and Wy, while the second set of terms will be determined solely by 



S = 0(A) 




p 
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amplitudes W^- Given the assumption that the leading order of A comes from W^, we have 
that the determinant of A will scale with jo as: 



Namely, the scaling of terms with diagonal elements in the determinant cannot be established 
without the detailed knowledge of its dependence on and df terms. However, the scaling 
of each off-diagonal component Akp{k) (where p{k) k, k = 1, . . . , R) will be determined 
only by the vertex amplitude Wy, and is dominant by assumption. As a consequence, the 
terms in det A which do not contain any diagonal elements are dominant and scale as O(A^), 
while the terms which do contain diagonal elements may scale with smaller power in jo and 
can be neglected in the limit jo — > oo. 

Once the scaling of det A has been established, we can employ some well-known results 
about the Schur complement matrix in order to establish the scaling of S. These results are 
summarized and proved in the form of Lemma 1 in Appendix [Dl 

Let the Hessian matrix A, its submatrix M and its Schur complement S scale as 



Note that M, being the submatrix of A, scales with jo with the same power —d as A. 
However, this cannot be assumed for the Schur complement S since there might be nontrivial 
cancellations between the leading terms in A and NM~^N'^ in f lTSj) . Consequently, the 
scaling power of 5* is —s. What we need to prove is that these cancellations do not happen, 
and that in fact s = d. 

Denote the ranks of M and S matrices as r and p, respectively. By part (b) of the Lemma 
1, we have 



Calculating the scaling order of the left-hand and right-hand sides, and using the usual 
properties of determinants, we easily see that 



det A = O(A^). 




det A = det M det S. 



1 



1 1 




which gives 



Rd = rd + ps. 



22 



By the part (a) of the Lemma, we have R = r + p. Using this to ehminate both R and r, 
the above equation reduces to 

p{s-d) = 0. 

Finally, by assumption of the theorem, matrix 5* is nonzero, which means that its rank p 
is nonzero. Therefore we conclude that s = d, which actually means that S = 0(A). As 
matrix 5* is a jj submatrix of S*, it scales in the same way as 5*. Consequently, 

S = 0(A) , 

which proves the theorem in the case when A is nondegenerate. 

If A has zero eigenvalues, the determinant equation above vanishes identically. However, 
in this case we can repeat the whole analysis in the same way, except that we need to use 
part (c) of the Lemma instead of part (b), bearing in mind that 0(54) = 1 (see Remark 
3 in Appendix [D]) . Namely, instead of analyzing the determinants of A, M and 5*, we can 
rotate the basis to represent these three matrices in the form 


















, M = 




, s = 




Ma 




M 




Ms 



and repeat the whole proof using the regular minors Ma, M and Mg instead. Note that as 
a consequence of the part (a) of the Lemma, the sum of dimensions of the zero-blocks of M 
and 5* must be equal to the dimension of the zero-block of A. These zero-blocks represent 
the kernel of A, and as discussed in the main text, appear as a consequence of continuous 
symmetries of the action (j5]). As was shown in section IVj they may safely be integrated out, 
and the Schur complement f|T8l) constructed from the regular part of M, i.e. the minors M 
and N. 

Again, since S* is a submatrix of S, if it is nonzero it scales with the same power as S, so 
consequently we have 

S = 0(A) , 

in the degenerate well. This completes the proof of the theorem. 

Appendix D: Properties of the Schur complement 



Here we establish some properties of the Schur complement that we have used in the 
proof of Theorem 1. These results can be found in 22|. However, one of the results, the 
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statement (c) below, is a new result, generalizing the statement (b). 

Lemma 1. Let A be a symmetric complex matrix of type n x n and let R be its rank. 
Let us split A into blocks as 

A N 

A = 

M 

where A is a J x J matrix, N is a J x r matrix, M is a r x r matrix and n = J + r. We 
will also assume that M is invertible, and that real and imaginary parts of A commute. 
Let us construct the Schur complement S (see which is a J x J matrix 

S = A- NM-^N^. 

Denote the rank of S as p. Then 

(a) R = r + p (Guttman rank additivity); 

(b) det A = (lets' detM (Schur determinant formula); 

(c) ifO<p<J, then 

det MA(det B^f = det M det Mg. (Dl) 



Here Ma one? are invertible Rx R and px p matrices, respectively. They are obtained 
by using orthogonal transformations which put A and S into a block- diagonal form 



A 



The B4 matrix will be explicitly constructed in the proof below. 













S = 




Ma 




Ms 



Proof. We start from the Aitken block diagonalization formula 22| and from now on we 



use / to denote a unit matrix of any size appropriate for its position in an equation: 



/ -NM-^ 
/ 



A N 
M 



I 





S 




M 



(D2) 



This equation can be verified by a direct multiplication of the left-hand side. Denoting the 
first matrix on the left as C, we can rewrite this identity in a compact form CAC^ = S(BM. 
The rank of the right-hand side is the sum of ranks of S and M, which amounts to p + r. 
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Since the rank of C is equal to its dimension n, the total rank of the product on the left-hand 
side is equal to the rank of A, so we easily obtain 

R = r + p, 

which completes the proof of part (a). 

Next, we take the determinant of (1D2I) . Since C is block-triangular, its determinant is a 
product of determinants of blocks on the diagonal, so we obtain det C = 1. The left-hand 
side is thus the product of determinants, det C det A det , and it is equal to det A because 
det = detC = 1. On the right-hand side we have a block-diagonal matrix, so that its 
determinant is equal to det S det M. Hence, 

det A = det^detM, 

which completes the proof of part (b). 

In order to prove (c), let O be a J x J orthogonal matrix which transforms S into a 
block-reduced form, 

OSO'^ = © M5 . 

Since p 7^ 0, matrix S has exactly p nonzero eigenvalues, which constitute Mg. Given that 
the eigenvalues of Mg are nonzero, it is invertible. The zero-block is of type z/ x z/, where 
V = J—p is the dimension of the null-space of S. By using O one can construct an orthogonal 
n X n matrix P = O ® I such that 

P {S ® P'^ = Q ® M-s® M. (D3) 

By using an analogous argument one can always construct an orthogonal n x n matrix 
such that 

Q^Ag = © Ma, 

which can be solved for A: 

A = g (0 © Ma) . (D4) 

The zero block comes from the null-space of A. It is of the size n — R, which is also equal 
to u, since n = J + r and R = r + p according to the part (a). 

Consider flD2|) . and multiply it by P from the left and by P^ from the right, and use 
f IDSP and (]D4|) to rewrite it in the form 



PCQ (0 © Ma) Q^C^P^ = 0®Ms®M. (D5) 
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Let us introduce the matrix B = PCQ and write it in the block form as 



B 



Bi B2 
B3 B4 



where the blocks Bi, B2, -B3 and i?4 are uxi/jUxR, Rxu and Rx R matrices, respectively. 
Substituting this into the left-hand side of (IDSp yields 



PCQ (0 © Ma) Q^C^P^ = B 




Ma 



B' 



B2MABJ B2MABJ 
B^MaB^ B^MaBJ 

By comparing (1D6P to the right-hand side of (IDSp . we obtain 



(D6) 



B2MABI B2MABJ 
B^MaB^ B^MaBJ 




M^ 
M 



(D7) 



Note that the zero-block of (lD7p is a. u x u matrix, which is also the B2MAB2 block. We 
then read off the following equations 



B^MaBI = Ms®M, 



B2MABI = 
B2MABI = 



(D8) 

(D9) 
(DIO) 



By taking the determinant of (IDSp . we finally obtain 

det MA(det B^f = det M det M^ . 

This establishes flDip and completes the proof of part (c). 

Given that M, Mg and Ma are all invertible, we have det B4 which means that B4 
is also invertible. By multiplying (]D9|) by (Bj)^^M^^ from the right, we obtain 



B2 = 0. 
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The equation (IDlOp now vanishes and does not provide any additional constraint. Therefore, 
the matrix B has the following form 



B = PCQ 



Bi 

B3 B4 



(Dll) 



End of proof. 

Remark 1. The A matrix from the main text has the form 

A N 
M 


which differs from the one in Lemma 1 by an additional zero-block. However, these additional 
zeroes are integrated out before the lemma is applied, and hence they do not affect any 
statements of lemma. 

Remark 2. The result (c) is a generalization of the result (b) to the case when A is a 
singular matrix. While the part (b) is in fact valid for singular matrices, it merely states 
that = and provides no information about nonsingular principal minors of A. The result 
(c) is more fine-grained, and provides precisely this nontrivial information about A. 

It was assumed in the part (c) that < p < J. \i p = J then A is a regular matrix, and 
hence the result (b) can be used. If p = 0, then 5 = 0, = J, and instead of f IDSp we obtain 



B^M^Bi = M 



and consequently 



det MA(det B^f = det M , 
In this case we can set P = I and obtain 



B = CQ 



Bi 

-B3 -B4 



for the matrix B. 



Remark 3. In Appendix O we use the results (b) and (c) to determine the leading 
jo-order of the Schur complement 5, knowing 0(A). However, it is necessary to show that 
-B4 is of order 0(1). In order to do this, note that 

det5 = detPdetOdetQ = ±1, 
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since P and Q are orthogonal matrices. On the other hand, from ( IDlip we know that 
det B = det Bi det i?4, so that we have 



det Bi det S4 = ±1 . (D12) 
Let us now assume that the blocks Bi and B4 are of order k and m in 1/jo, respectively 

k,m>0, D,Er^O{l). 

The numbers k and m cannot be negative since the whole B matrix must be of order 0(1). 
Namely, the matrices P and Q are orthogonal, and consequently all their elements are 
bounded above by 1. Thus P and Q are 0(1). The matrix C is also 0(1), since A and 
consequently M, A^, are all of the same order. Therefore, B = PCQ ~ 0(1). 
Since i?i is a x matrix and B4 is a R x R matrix, we have 

detSi = -^detD + o(-^) , (D13a) 

det B, = -^detE + (-^) . (D13b) 
Jo \Jo J 

Substituting (1D13P back into (]D12p we obtain the consistency equation 

kv + mR = . 

Since both i/, i? > while k,m > 0, the only solution of this equation is k = m = 0. 
Therefore 

det 54 - ^4 ~ 0(1) . 

In the case when z/ = the A matrix is regular and instead of the part (c) we use the 
part (b) of Lemma 1. However, the part (b) does not involve det-B4, so that we need the 
above result only for z/ > 0. 

Appendix E: Vertex amplitude contribution to the Hessian matrix 

Here we give the explicit formulae for the terms on the right-hand side of ( !25l) . These 
terms are calculated by directly substituting the vertex asymptotics ( ]B1|) into equation ( 12^ 
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and differentiating. It is important to note that in the expression (IBip the x-dependence is 
in the coefficients and A^_, as well as in the Regge action Sr. However, the scaling of 
N± is different than that of Sr. The former scale as 0(1) while the latter scales as 0(jo) in 
the limit jo ~^ oo- 

The 0(1) term in has already been quoted in the text in equation and we 

repeat it here for completeness: 



A 



(0) 
vab 



The 0(jQ ^) term is given as: 



vab 



+2. 



jia6a,b£vQ 



(v) 
b 



The 0(jQ ^) term is given as: 



A 



(2) 

vab 



5, 



a,b(iv 



Appendix F: Comparison with the Euclidean results 

It is instructive to compare our results with the asymptotic analysis of the Euclidean 
ELPR/FK state-sum kernel given in ISj. There are several differences in the setting between 
the approach of ISj and the one taken in this paper. Since they do not consider the boundary 
wavefunction, they do not have the insertion functions. Next, they integrate the state sum 
over all variables except the spins, and obtain 

where WA(j/) is the stat-sum kernel, A is the triangulation dual to a and the face amplitude 
is a product of two df terms. 
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One of the main results of 18| is the asymptotic large-spin expression for the nondegen- 
erate part of this kernel 

where is the large parameter, ca is the function of spins, but not of and Sr is 



see equations (87) and (101) in 18j. Note that 0/ is the sum of all dihedral angles around 
a face /. The action Sr is constructed for the triangulation A, with non-scaled spins, and 
the large parameter is written explicitly in front of it in W/\. In our notation, the above 
expression can be rewritten as 

W^'^{j) = A{j)cos{Sr), j^oo, 

where now Sr is constructed with the scaled spins j, while the amplitude is denoted simply 
by A{j). 

Despite all the differences in the two setups, there is a rather simple generic relation be- 
tween the kernel W^a(j) and our boundary wavefunction ([2]). The quantity that corresponds 
to ([2]) can be constructed from the kernel W^{j) in the following way. First we choose the 
triangulation A so that it has a boundary. The dual of A will be the 2-complex a, while 
the dual of the boundary will be a 1-complex 7 = da. Next, we split the face labels jf 
into the boundary and internal labels, and denote them jf and kf, respectively. Then the 
wavefunction is given as 

kf f 

Here we have not introduced the insertion functions fii on the boundary, and the face am- 
phtude is quadratic in spins. 

In the limit where both internal and boundary spins are large, one can approximate the 
sum with an integral and write the wavefunction as 

= [dk llid+d_)Aij, k) cos(5^(j, k)) . (Fl) 

In order to evaluate the integral over k in the large- spin limit, we would like to approximate 
the cosine with a Gaussian in the neighborhood of each of its stationary points, and employ 
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the stationary-point method. This technique was used for the asymptotic analysis of the 



wavefunction for the Euchdean canonical LQG jo], which is given by a similar state sum 
as (IFip . The cosine has infinitely many stationary points, but we can assume that the 
amplitude A{j, k) is peaked around only one of them (otherwise the asymptotics will never 
be a Gaussian function; also one could introduce the insertion functions which would single 
out one stationary point). Denote this stationary point as (j*, k*). In the neighborhood of 
this point, the cosine can be approximated with a Gaussian via the formula 

. Using the fact that x* is the stationary 



where x = {j,k) and A^^ = log cos Sr{x) 

point of the cosine and the fact that the variation of Sr w.r.t. the angles 0/ is identically 
zero, an explicit evaluation of Aab gives 

Aab = -Qa&b = 0(1) . 

At this point one can split the matrix A into jj, jk and kk blocks, like in Eq. ( !T6|l 
and perform the integration over fc-spins. The result will be a Gaussian over the remaining 
j-spins, 

where S is the Schur complement of A. Since the order of A is 0(1), so will be the order of 
the matrix S, confirming our general result ( p7|) . 



[1] C. Rovelli, Quantum Gravity, Cambridge University Press, Cambridge (2004). 

[2] A. Mikovic, Class. Quant. Grav. 21 3909 (2004), |arXiv : gr-qc/040402fj . Errata: Class. 

Quant. Grav. 23 5459 (2006), a rXiv : gr-qc/ 060608lt 
[3] C. Rovelli, Phys. Rev. Lett. 97 151301 (2006), a rXiv:gr-qc/0508124 

[4] E. Bianchi, L. Modesto, C. Rovelli and S, Speziale, Class. Quant. Grav. 23 6989 (2006), 

|arXiv : gr-qc/0604044y 
[5] E. Alesci and C. Rovelli, Phys. Rev. D 77 044024 (2008). larTiv : 071 1 . 1281 
[6] A. Mikovic, Fortschr. Phys. 56 475 (2008) . larXiv : 0706 . 04661 

[7] E. Bianchi, E. Magliaro and C. Perini, Nucl. Phys. B822 245 (2009) . [arXiv : 0905 . 40521 



31 



B. Dittrich, L. Freidel and S. Speziale, Phys. Rev. D 76 104020 f 20071 . liiiTiv : 0707 . 45 131 
A. Mikovic and M. Vojinovic, to appear in Adv. Theor. Math. Phvs.. [arXiv: 1005. 18661 

J. Engle, E. Livine, R. Pereira and C. Rovelli, Nud. Phys. B799 136 (2008), 
larXiv: 0711 .01461 

L. Freidel and K. Krasnov, Class. Quant. Grav. 25 125018 (2008) . liEoCiv : 0708 . 15951 

J. B. Hartle and S. W. Hawking, Phys. Rev. D 28, 2960 (1983). 

E. R. Livine and S. Speziale, Phys. Rev. D 76 084028 f 20071 . E?Xiv : 0705 . 06711 

C. Rovelli, l arXiv: 1004. 1780 1 

E. Bianchi, D. Regoli and C. Rovelli, larXiv: 1005.07641 

J. Barrett, R. Dowdall, W. Fairbairn, H. Gomes and F. Hellmann, J. Math. Phys. 50 112504 
(2009), larXiv: 0902. 11701 

L. Hormander, The Analysis of Linear Partial Dierential Operators I, Spring er-Verlag, Berlin 
(1983). 

F. Conrady and L. Freidel, Phys. Rev. D 78 104023 (2008), arXiv: 0809. 22801 

E. Magliaro and C. Perini, larXiv: 1105.02161 

J. Barrett, W. Fairbairn and F. Hellmann, larXiv: 0912. 49071 

J. Barrett, R. Dowdall, W. Fairbairn, F. Hellmann and R. Pereira, larXiv: 0907. 2440 | 

F. Zhang, The Schur Complement and its Applications, Springer, New York (2005). 

J. Baez, J. Christensen and G. Egan, Class. Quant. Crav. 19 6489 (2002), 
|arXiv:gr- qc/02080J^ 



J. Barrett and C. Steele, Class. Quant. Crav. 20 1341 (2003), [arXiv:gr-qc/0209023 

J. Barrett, R. Dowdall, W. Fairbairn, H. Gomes, F. Hellmann and R. Pereira, 

larXiv: 1003. 18861 



32 



